TODA HIERARCHY, HURWITZ NUMBERS 
AND CONFORMAL DYNAMICS 



S.M. NATANZON AND A.V. ZABRODIN 



Abstract. We explicitly construct the series expansion for a certain class of solutions to 
the 2D Toda hierarchy in the zero dispersion limit, which we call symmetric solutions. We 
express the Taylor coefficients through some universal combinatorial constants and find 
recurrence relations for them. These results are used to obtain new formulas for the genus 
double Hurwitz numbers. They can also serve as a starting point for a constructive 
approach to the Riemann mapping problem and the inverse potential problem in 2D. 
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1. Introduction 

The dispersionless 2D Toda lattice (2DTL) hierarchy was introduced by Takasaki and 
Takebe |28[ 129] . It can be represented in two equivalent ways: in the Lax-Sato form or 
in the Hirota form. In this paper we use the latter formulation. The dispersionless 2DTL 
hierarchy is an infinite system of differential equations 

(1) {z - Oexp{D{z)D{C)F) = zexp{-doD{z)F) - ^exp{-doD{OF), 

(2) {z-Oexp{D{z)D{OF) = zexp{-doD{z)F) - ^exp{-doD{OF), 
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(3) 



1 - exp{-Diz)D{OF) = lexp(9o(9o + D{z) + 5(0)^) 



for the function F = F{to,t,t), where t = {ti,t2, ■ ■ ■}, t = {ti,t2, . . .} are two infinite 
sets of time variables and 

(4) D{z) = J2'-^d,, D{z) = Y,'-^dk. 

k>l k>l 

Hereafter we abbreviate dk = d/dtk, dk = d/dik- Differential equations of the hierarchy 
are obtained by expanding equations ([I])-® in powers of z, ^. These equations are known 
to be connected with different branches of mathematics and mathematical physics. 

Solutions to equations ([I])-© such that dkF\tQ = dkF\tg = for all A; > 1 form 
an especially important class. (By gltoito) we denote the restriction of a function 
g{to,ti,ti,t2,i2, . . .) to the line ti = ii = t2 = h = ■ ■ ■ = 0.) We call them symmet- 
ric solutions. This class contains, in particular, the c = 1 string solution [21 El El ED]- 

In Section 2 we describe all formal symmetric solutions of the dispersionless 2DTL hier- 
archy in the form of a Taylor series. We prove that they are fully defined by the restriction 
F\tf^ (a function of one variable). The Taylor coefficients are expressed through some uni- 
versal combinatorial constants Nf^^^^^-j / '^i • ■ ■ j which depend on two Young diagrams 

y 1 • • • J 

A, A and two sequences of natural numbers {sj}, {vi}. Moreover, we find recurrence for- 
mulas for A/'i^AlA) I "^^ ■ ■ ■ "^"^ I . Our method is an extension of the method developed in 

[22| |23] for the formal c = 1 string solution. Convergence of the Taylor series for this case 
was investigated in [llj. 

In Section 3 we apply these results to the double Hurwitz numbers for coverings of 
genus 0. This application is based on the fact that the generating function for connected 
double Hurwitz numbers of genus is a symmetric solution of the dispersionless 2DTL 
hierarchy, which follows from Okounkov's result [25] and its dispersionless limit discussed 
by Takasaki [27J. 

The connected genus double Hurwitz number Hq{A\A) is the number of non- 
equivalent rational functions with pre-images of and oo of fixed topological types given 
by Young diagrams A = [ki, . . . , ki] and A = [ki, . . . , kg] (with the condition on degrees 
|A| = |A|) and some fixed simple critical values (finite and non-zero). For rational func- 
tions of degree d with simple finite critical values this number was found by Hurwitz 



d-2 



where Yli^i ~ ^ ^'^'^ ^(^) order of the automorphism group of rows for the Young 
diagram A. 

Further results are based on the moduli space methods [TB]. In [T7] it was found, in 
particular, that 

i/o(N|A) = i!^n--. 
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Some formulas for Ho{[ni,n2]\A) and Ho{[ni,n2,ns\\A) were found in [HI [26]. It was also 
proved that the Hurwitz numbers iJo([[^i; • • • ? ^n]]|[^i; ■ ■ ■ , ^n]) depend piecewise polyno- 
mially on ki, kn,ki, ... ,kn [6J. The domains of polynomiality (chambers) were char- 
acterized and relations between polynomials in neighboring chambers were found [26j . 

In Section 3, by comparing the generating function for the genus double Hurwitz num- 
bers with the formulas for symmetric solutions, we give a representation of the Hurwitz 

numbers through the combinatorial coefficients N(^^^^^ ' ' These formulas ap- 

pear to be new and useful. In particular, they make explicit some general properties of the 
genus double Hurwitz numbers (for example, the piecewise polynomiality). Moreover, 
they are well-suited for direct calculations with the help of the computer and allow one 
to obtain closed formulas for ifo([^i; • • • ? ^n] I [ki, . . . ,kn]) for any n, n. We illustrate this 
method by two examples: Ho{[n]\[ki, . . . , k^]) and Hodki, k2]\[ki, /C2]), where the calcula- 
tions can be done by hands. Other consequences of our general formulas will be discussed 
elsewhere. 

Another application, addressed in Section 4, is the class of problems which we call 
conformal dynamics. In papers [T^ [T^ ITS] it was shown that some problems of 
complex analysis in 2D, such as conformal mapping, Dirichlet boundary value problem 
and 2D inverse potential problem have a hidden integrable structure, which, for simply- 
connected domains, is the dispersionless 2DTL hierarchy. In a more general case it is the 
universal Whitham hierarchy introduced by Krichever in [T3 | [T ^ . The hierarchical times 
tk are suitably defined harmonic moments of the domain and their complex conjugates, 
with to being proportional to the area of the domain. One can construct the dispersionless 
tau-function F = F{tQ, {tk}, {ik}) which contains all the information about the conformal 
bijection of any domain with given moments to the unit disk. The function F satisfies the 
dispersionless version of the Hirota equations for the 2DTL hierarchy given by ([I])- ([3]). 

Further, in [35] it was argued that any non-degenerate solution of the hierarchy, with 
certain reality conditions imposed, can be given a similar geometric meaning. Such so- 
lutions are parameterized by a function cr{z, z) of two variables which has the meaning 
of a density, or conformal metric, in the complex plane. The moments should be now 
defined as integrals of powers of z with this density. The integral representation for the 
dispersionless tau-function also changes accordingly but the formulas which express the 
conformal map through its second order derivatives do not depend on a. In other words, 
the Toda dynamics encodes the shape dependence of the conformal mapping, which we 
call the conformal dynamics. In the context of the conformal dynamics, the symmetric 
solutions correspond to the axially symmetric functions a (i.e., depending only on \z\'^). 

Some important examples of conformal dynamics on symmetric background are con- 
sidered in Section 5. Using the approach of Section 2, one can write explicit formulas 
for the corresponding symmetric solutions of the dispersionless 2DTL hierarchy. The 
corresponding conformal dynamics provides an effectivization of the Riemann mapping 
theorem. Namely, for any domain characterized by its moments, our formulas allow one 
to find the conformal map from this domain to the unit disk with any given precision. 
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2. Formal symmetric solutions of dispersionless 2D Toda hierarchy 



Let A = [yUi, /i2, . . . , /i£] be the Young diagram with ^ = i{A) rows of non-zero lengths 
/Ui > At2 > . . . > /i£ > 0, and similarly for A = [fxi, fi2, ■ ■ ■ yfie], with i = i{A). We identify 
A with the partition of the number |A| := /ii + . . . + /i^ into the i non-zero parts /Xj. 
Another convenient notation is A = (^l™-i2™2 _ _ f"^r _ _ which means that exactly 
parts of the partition A have length i: rrii = card{j : = i}. Put a(A) = mi! . . .m^! 
and p(A) = /ii . . .//^(a). 

Given the two sets of time variables, t = {ti,t2, . . .}, t = {ii,t2,---} as before, set 
tA = tf_i^ . . . t^i,, tA = ^/ii • • • ifif For the empty diagram we put = = 1. We say that 
a Taylor series of the form 

(5) Fito, t, t) = ^ F(A| A|to) tA^A, 

A, A 

where the sum is over all pairs of Young diagrams including empty ones, is a formal 
solution to the dispersionless 2DTL hierarchy if its substitution to equations ([I])-© gives 
the identity. In a more explicit but less compact notation the series reads 

F(to, t, t) = F(to) + ^ F(/Xi,/i2, . . . ,/i£|/il,/i2, • • • ,/i£|i!:o)i^Mii^M2 • • • '^Mfc^/il^/i2 • • -^M^-- 

M1>P2>--->M|" 

Here F{to) = F(0|0|to). 

Remark. In this section the bar does not mean the complex conjugation and tk, ik are 
regarded as arbitrary formal variables. 

2.1. Taylor expansion of symmetric solutions. Let gltoito) denote the restriction of 
a function g{tQ, ti, ti, t2, ^2, • • • ) to the line ti = ti = ^2 = ^2 = • • • = 0. Formal solutions 
such that dkF\tQ = dkF\tQ = for all /c > 1 will be called symmetric. 

Our first goal is to prove that symmetric formal solutions are fully determined by the 
function F|j(,(to) := F{tQ) which can be an arbitrary twice differentiable function of tg- 
Moreover, we prove that $(A|A|to) is a differential polynomial in /(to) = 6xp(F|5.'p(to)) 
(i.e., a polynomial in f{to),f'{to),f"{tQ),...) with universal coefficients which will be 
found. For a particular function /(to) they were found in [22j. In the general case the 
arguments are similar. We split the proof into several lemmas. 

Lemma-Definition 2.1. (|22). Lemma 3.3) For any formal solution F the following 
relation holds 

= — ~ Tij{pi, ...,pm) dodp,F. . . dodp^F, 

m>0 Pl---Pm 

PlH t-pm.—i+i 

where 

TijiPl, ■■■,Pm)= -j^i — \ Yp'' Y ■ ■ ■ ^ Y p' 

k>o,ni>0 1- ■ ■ ■ k- \ i=l i=ni + l i=ni+...+rifc_i + l 

h'T'/j.— TTl ^ 

and Pij{ri, . . . ,rm) is the number of sequences of positive integers (ii, . . . ,im), {ji, ■ ■ ■ ,jm) 

such that ii + . . . + im = i, ji + . . . + jm = j and rk = ik+ jk- 
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An analog of this lemma for the dispersionless KP hierarchy can be found in |24) . Here 
and below all indices and variables with indices like i,p, s, i, n, r, a, b with or without bar 
are positive integers. 

By induction one can prove 

Lemma-Definition 2.2. ([22], Lemma 3.4-) For any formal solution F the following 
relation holds 

oo ■ ■ / \ 

m=l n + ' ■+=m=n+...+«fc 1 • • • m V 1 m/ 

ei + ---+em=m + k-2 

where 



-^ll...lfc 



Si ... s. 



Si . . . 
^1 • • • 



a, -IV.. 



otherwise, 



^i-.-C; ^ -1)! 



m / "^1 • • • Si — I S Sj-^-l . . . Sm \ rp j \ 

with 

s = Si + Si+i + ■ ■ ■ + - ifc > 0, £ = (^i - 1) + ■ ■ ■ + - 1) > 0. 

From now on we consider only symmetric formal solutions F . Put / := exp(F|"|^^ 
Lemma 2.1. If F is a symmetric formal solution, then 

fori^ j , 



Proof. From 9o9,fcF|^^= (9o9fcF|^^= for A; > it follows that 

exp{do{do + /^(^) + ^(O)i^) = exp(9o^F|,J = exp(F|;'J = /. 



to 



Moreover, from 1 - e^^^^)^®-^ = ;2-i^-ie9o(^o+i?(2)+5(0)F ^^^^ 

oo ^ 

-D(z)D(OF|^^= log (1 - z-'t'f) = - E ^^"'rV'. 

fc=i 

Therefore, 9i9jF|^^= for i 7^ j and didiF\^^= if^. □ 
Lemma 2.2. The following relations hold 

if ii + . . . + ik i 

ii . . . ikdo'\f') if ii + ... + ik = i. 



= didi^ . . . di^F 

to 



to 



Proof. The differentials d and d enter the Toda equations in a symmetric way. This gives 
the first equality. Moreover, according to Lemma-Definition \2.2\ we have 

didi,di2 ■ ■ ■ \F = . ^^•••^^ do'^didi^+,„+i^F 
ti + . . . + tk 
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oo . . . ^ 

■^"^ OT . . . bm \ ^-1 • • • '-m / 



m=2 si + ... + Sm=ij^+. .. + !(, 
^l + ...^m="i+fe-2 



This equality and Lemma [2TT] gives the second equahty in the assertion of Lemma I2l2l □ 

Given a sequence of natural numbers {ai, . . . , a„}, we say that its representation as a 
union of non-intersecting non-empty subsequences {&{, 6^ ,} such that 5{ + . . . + 6^ . = Sj, 

{ai, . . . , a„} = U • • • > is a partition of {oi, . . . , a„} o/ i/ie t?/]9e ( ^ ' ' ' • 

j=l \'>'l ■ ■ ■ ll'mj 

m 

Clearly, the sequence (ni , . . . partition of v in the usual sense: ; 



n,: = t;. 

i=l 



Lemma-Definition 2.3. The following relation holds: 



to _ _ , . . 

m=l siH-... + sm=ii + ... + sj.=ii + ... + ij \tl 
r^H hrm— Tn+fc + fc — 2 



where 

,7 / si . . . \ _ ti • • • ffc n - ■ ■ ik sr^ rj. 



n ...jfc 



ri...rm)~ Si...Sm ^ Ui -m + 1 .. .r^-n^ + 1 



and the summation is carried over all partitions of the set {zi, ife} of type 



S\ . . . Sff 

ni . . .n„ 



Proof. According to Lemma-Definition \2.2\ 



di,...d,^\...di F 



00 . . / s 



m=l S]^ + ... + Sm=il + . 

Zl + ... + im— + — 2 



n • • • ^A: ^ / Si . . . S 



E E E 7^^'-..-M/, / 



m=l si + ... + Sm=il+-.- + ij; 

^X + .-.+Zm — m + fc — 2 



where the interior summation is carried over all partitions of of type 

S\ . . . 



. Therefore, according to Lemma [2 ■2[ 



di, . . . d.^di^ . . . d-, F 



m = 1 Si +...4-s^=i> +... + 17. '^1 • • • '^m, V \ ^ / / 



m=l si + ... + Sm=il + --- + ik 
l-l + ...+lm=rn + k — 2 



E E 



St . . . Sr 

m=l s]^+...+sm=ii+-.-+»j. 

rj + .-.+rm— m+fc + fc — 2 
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□ 



where the interior summation is carried over all partitions of the set {ii, . . . ,2^} of type 
Si . . . 
rii. . .rim 

Lemma-Definition 12.31 implies 

Theorem 2.1. Any symmetric formal solution to the dispersionless 2DTL has the form 

(6) F = F{to) + Yl E ^(A|A) 5o"(r)---55'"(r'")'^AtA, 

A| = |A| si + ...+sm = |A| _ \ ■ ■ ■ m/ 

ri + ... + rm=m+e(A) + l{A)-2 

where /(to) = exp{d^F{to)) , 

and . . . , /i^] = A, [fii, ...,fxt,] = A. 

Theorem 12.11 implies that the symmetric solutions are parameterized by a function of 
one variable. The Taylor coefficients in the right hand side of equation are differential 

Si . . .Sr 

^) ■ 







... 


■ • \ 




■■^^k\ 




■ ■ '^m J 




■■P-kJ 







polynomials in the function /(to) = e"^"^*"^ with universal coefficients A/^^^i^^ 



Ti . . .r. 



m 



S\ . . . Sm 

Explicit formulas for the simplest coefficients A''(a|a) ( ^ ' ] given below. 



ri . . . r 



m 



2.2. Explicit formulas for some coefficients. 

Lemma 2.3. Let A, A be two Young diagrams such that \A\ = \A\ = d. Then 

\A\ \ p{A)p{A) 



^{A|A) ^£(A) +£(A) -2)- da{A)a{A) 



and Ai"(^|A) = otherwise. 

Proof. According to our definitions Tjj ^ '^'^^ = + j) = + j) = 1 



Tij { ] = in other cases. Therefore, 



J-ii...ik\ k — 1 ) ~ ^«i-«fc-i 1 k — 2 I + ■ • • + — J- 

and Tj^ = in other cases. Let the sequences {ii, ■ ■ ■ ,ie(A)), iH,---,h{A)) 

be any permutations of [yUi, . . . , /i^], [fii, . . . , fte] respectively, where [fii, . . . , fie] = A, 
[/2i, ...,fii] = A. Then 

N . . ( I^L \ _ p(A)p(A) ^ / |A| ^_p(A)p(A) 



\£(A) +£(A) -2; |A| I £(A) - 1 / |A| 
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11 . . . 



and A^/. \ ( 1 = in other cases. Hence 



^V(A|A) ^£(A)+£(A)-2y' - a(A)a(A) ^V!--^-) 1^£(A) + £( A) - 2^1 " |A|a(A)a(A) 



*i . . . 

and A''(A|A) ("^1=0 otherwise. □ 



Lemma 2.4. 



%) (n) = %) (u) = - 2.(|....!fM|,...j) -"'°fe-'-'-'''> 



n«2 



and N ( "^""^"^^ 1 = otherwise. 



rir2 



Proof. The number Pij{ri,r2) is the number of solutions to the system of equation 
ii + h = h ji+j2 = j, H+ji = ri, 12+32 = r2 in positive integer numbers ii, ^2, ji, J2- This 
system is equivalent to the system ii = ri — ji, i2 = i+ji—^i, J2 = j — ji, i+ j = ri + r2. 
Hence Pij{ri, = if z + j 7^ ri + r2. 

In other cases the number of solutions is the number of positive ji such that ri —ji > 0, 
i + ji - n > 0, j - ji > or ji < ri, ji < j, ji > ri - i. By symmetry in pairs i,j 
and ri , r2 it is enough to consider the case j < i and ri < r2 . In this case the number of 
solutions is min{ri,j} — 1. Thus Py(ri,r2) = min{i, j, ri, r2} — 1. 

By our definition 

Tiui2 (^ll^^ =Tij{si,S2) = -^(^Pij{si + S2) + ^Pij{si,S2)^ = -^min{i, j,si,S2} 
if i + j = si + S2 and Ti-^^i^ f ^^f^ ) = in other cases. 



By Lemma-Definition 12.31 we have 

N fsiS2\ ^iji2Mi2Sr^rj, f Si S2 

ff^Q\rir2j S1S2 ^ ^i^^Vri-m + l r2-n2 + l 

where the summation is carried over all partitions of the set {ii, ^2} of type {^^^ ) • Such 
partitions exist only if either %x = Si, 12 = S2 or ii = S2, 12 = Si. Therefore, 

nr lili2\ nr A2«l\ «1^2^1^2^ Al ^2^ ^1^2 . r. . ^ 

^^i^lA [ 11 J ^^(i^iA \ 11 J " 1 ) = -— minOi,^2,^i,^2}) 



and in other cases. Finally we obtain 

N\^.^2^M ^ 11 J - N^nMM ^ 11 J - -2a([.i,Z2])a(pij2])^'^'''^''''"'''''^^ 
and ATf^ii^llPJa] Wirl) ^ ° o'therwise. □ 
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3. Double Hurwitz numbers 



3.1. Generating function for double Hurwitz numbers of genus 0. Given a mero- 
morphic function : i7 — )■ C of degree d on a connected Riemann surface fl, one can 
associate with each point p G C = C U {oo} a Young diagram A = A{ip,p) = [fii, . . . , fii] 
such that |A| = deg ip = d and /ij equals the degree of the map if at the point of 
the complete pre-image ^~^{p) = {p^,...,p^}. Informally speaking, /ij is the number 
of sheets of the covering that are glued together at the point lying above p. We say 
that meromorphic functions (f : Q ^ C and : f2' — )■ C are equivalent if there exists a 
biholomorphic map f : Q Q' such that (p = f'f'. 

Fix different points pi, . . . ,pk and Young diagrams Ai, . . . , A^ such that |Aj| = d for 
alH = 1, . . . , fc. Consider the sum 

</peCovd(n,{Ai,...,Afc}) ' 

where |Aut(y9)| is the order of the automorphism group the covering ip, 
CoVrf(f2, {Ai, . . . , Afc}) is the set of biholomorphic equivalence classes of meromorphic 
functions (p : Q ^ C such that A{(p,pi) = Aj and A{(p,p) = (l'^) at the other points. 
This sum does not depend on positions of the points pi. It is called the (connected) 
Hurwitz number. See [16] for a review. 



The Hurwitz numbers 

H,^i{A\A) = H{A, A, {l''-'2'), . . . , (l'^-22i)) 



I 



for functions of degree d with arbitrary ramification types A, A at two fixed points and 
I simple critical values of the type (l'^~^2^) are called double Hurwitz numbers. In fact 
the notation Hd^i[A\A) is redundant because this number is zero unless d = \A\ = \A\. 
The genus g oi Q and the number / are connected by the Riemann-Hurwitz relation 
2g — 2 = I — £{A) — i{A). Introduce two infinite system of variables t = {ti,t2, • ■ •}, 
t = {^1,^2; • • as in the previous section, and consider the generating function 

nl ^(A) ^(A) 

(8) ^f3,Q,t,t) = J2jJ2Q' E H,,l{A,A)l[^ia,.llfi^i,^, 

l>0 ' d>l |A| = |A|=d «=1 «=1 

where A = [/xi, . . . , //^(a)], A = [fn, . . . , fie^A)]- 

In order to extract the contribution of genus g surfaces the following trick can be used 
(see, e.g., [T] for the case of ordinary Hurwitz numbers). Let us rescale {tfc}? {ik} and /3 by 
introducing a new parameter h as tj. ^ t^/ h, f3 ^ hfS and consider the modified generating 
function /3, Q, t, t) := h?^{hf3,Q,t/h,t/h), then the series ([8]) having regard to the 
Riemann-Hurwitz formula acquires the form of the topological expansion 

(9) (3, Q, t,t) = J2 Q, t, t), 
where 



(10) "^^ = ^1 fflA)+£(A)+2(7-2)! ^'^'^^^>^(^>"^^^^^'^^n^-^A^-nM 



^(A) ^(A) 

r, 

Mi 



counts the connected coverings of genus g. In particular 

/32(£(A)+£(A)-2 



d>l |A| = |A|=d 

is the generating function for the numbers of the ramified coverings C — 
In [25] Okounkov has proved that the (dispersionfull) tau-function 

r„(t,t) = eii''"("+i)(2"+i)Q5"("+i)exp(^$(/3,e^("+5)Q,t,t) 

solves the 2DTL hierarchy of Ueno and Takasaki [32]. The dispersionless limit was char- 
acterized in terms of string equations by Takasaki [27]. The procedure of passing to the 
dispersionless limit [29j is equivalent to extracting the contribution of genus surfaces. 
Thus the function 
(12) 

F{/3, Q, to, t, t) = ^ + ^ log Q + $o(/3, Qe^*« , t, t) 

d>l |A| = |A|=d ^" V V / V / / j=i j=i 

satisfies the dispersionless 2DTL hierarchy from the previous section. Moreover, from the 
structure of the series we see that it is a symmetric solution. 

3.2. Formulas for double Hurwitz numbers. Now we can use Theorem 12.11 to cal- 
culate the genus double Hurwitz numbers. Note that at genus the Hurwitz num- 
ber Hdj{A\A) depends on the diagrams A, A only (the numbers d,l are restored as 
d = \ A\ = |A|, / = i{A) + i{A) - 2). We will write Ho{A\A) = Hd,i{A\A). 

Theorem 3.1. The genus double Hurwitz number Hq{A\A) is given by 

where the sum is carried over all matrices I ' ' ' \ such that s-i + . . . +Sm = \A\ and 

\ri.--rmj 

ri + . . . +r^ = m+£(A)+£(A)-2 
Proof. The function 

^3 JA\to H (A A) '^^^ '^^^ 

F{wo,t,t)=j+ E (£(A)+/(A)-2)! n(/^-^^^)n(/^-^'Aj 



|A| = |A 



i=l i=l 



is a symmetric solution of the dispersionless 2DTL hierarchy with /(to) = Therefore, 
according to Theorem 12. 1^ 

^(i,i,to,t,t) = |+ j2 E ^(AiA)(;;:::';;)5o"(e^^*")---5S'"(e^'"*'')'^AtA 

IAI_IAI ., ,. lAI ^ / 



|A| = |A| si + ... + s„ = |A| 

+ . . .+r™=m+^(A)+£(A)-2 
10 



Comparing this with the previous formula, we find that 

Sl + . . . + Sm = |A 

r-1 + . . . + r„ = m + £(A) + ^(A) - 2 

which is the assertion of the theorem. □ 
Remark. Theorem 13.11 and the geometric definition of Hurwitz numbers give some 
nontrivial combinatorial relations. First, the sum ^ s'^ . . . s''^ ^(a|A) | "^^ ' ' ' "^"^ J must 

be positive despite the fact that the numbers A/^^^i^^ l^i---Sm\ have different signs. 
Second, the equality i/o(^|^) = -f^o(^l^) implies the identity 

which is non-trivial because all the definitions depend on the order of the Young diagrams. 

It follows from Theorem 13.11 that the formulas for A/^j-^i^) ' ' give explicit ex- 

pressions for the genus double Hurwitz numbers. Let us consider some examples. 

Corollary 3.1. The number of polynomials of degree n with a single critical value at 
of the type A is 

J =0 for m > 1. Hence, 

Til . . . f^m J 

according to Theorem 13.11 and Lemma 12. 3[ 



□ 



H fAlf.n (^(^) - |AK(A)-1 P(A)P(M) fflA) - 1)! .(A)-2 



Corollary 3.2. Let ii + 12 = ii + 12 = d. Then 

TT (I- -iir ^ c?- min{il,^2,^l,^2} 

^^o([h,^2] 1^1,^2]) = 2 — ^ 

Proof. Put A = [^1,^2] 7 A = [«i,«2]- It follows from our definition that 
^(A|A) I ■ ■ ■ ^™ j = for m > 2. According to Theorem 13.11 and Lemmas 12.31 12.41 
we thus have: 

(f(A)+£(A)-2)! ^ 
= p(A)p(A) l^'^---'- ^r, . . . J 

(£(A)+£(A)-2)! / _ / |A| V 277 TV ('^A 



p(A)p(A) 
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(£(A)+£(A)-2)! /|^|,(A)+nA)-2 P(A)p(A) ^l^^li2 ... ■ ^ ^ . \ 

P(A)MA) V' ' |A|a(A)a(A) a(A)a(A) ^ i' i' 

(^(A) +£(A)-2)! / 3 \ |A| - min{ii,i2,«i,^2} 



a(A)a(A) V ' ^ ^' ^' ^' c^(A)(7(A) 

□ 



4. CONFORMAL DYNAMICS 

The result of section 2 takes on a geometric significance when the Toda times tfc, are 
identified with (complex conjugate) moments of simply-connected domains in the complex 
plane with smooth boundary. In this case the dispersionless 2DTL dynamics encodes the 
shape dependence of the conformal mapping of such a domain to some fixed reference 
domain, thus providing an effectivization of the Riemann mapping theorem. We call it 
the conformal dynamics. 

4.1. Local coordinates in the space of simply-connected domains. Let D C C be 

a compact simply-connected domain whose boundary is a smooth curve 7 = dD and let 
D'^ = C \ D be its complement in the Riemann sphere C. Without loss of generality, it is 
convenient to assume that D 9 (and so D"^ 3 00). 

Fix a real-analytic and real- valued function f/(-2, ^) in C* = C \ {0} such that 

a{z,z) := ddU{z,z) > 0. 
We introduce the set of moments of the domain D*^ as follows: 

tk = ^—r <P z~''dU(z,z) dz = - [ [ z'^aiz,z)d^z , /c > 1, 
27r«fc Tck J Joe 

where d'^z = dxdy. In general they are complex numbers. We claim that together with 
the real parameter 

to = ^— (p dU {z, z) dz = — 11 (7{z, z) d'^z 

they form a set of local coordinates in the space "H of such domains D (or D'^). The 
function cr plays the role of a background charge density in the complex plane. 

First we prove the local uniqueness of a domain with given moments [4', T!S]. 

Proposition 4.1. Any one-parameter deformation D(t) of D = D(0) with some real 
parameter t such that all tk are preserved, dttk = 0, k > 0, is trivial. 

Proof. The proof is a modification of the one presented in [15j for the case a{z,z) = 1. It 
is based on the following two facts. 

• The difference of the boundary values dtC^{C,)dC, of the t-derivative of the Cauchy 
integral 

dz j; Unicode 



c{z)dz = — f \ UdCO = d^UiCC), 

^71% J ^ (, Z 
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is a purely imaginary diffential on 7. Indeed, let (^{9, t) be a parameterization of 
the curve ^{t). Taking the ^-derivative of the Cauchy integral and integrating by 
parts, one gets 

d,C{z)dz - ^-^ (C^J^^ 

dz / CtCe - CtCe 



2m 



f^-^^^-^a{ae,t),ae,t))de 



Hence, (9iC+(C)-9,C-(C))dC = <7{C,C){dtCdC-dtCdO = 2iaIm(dtCdC) is indeed 
purely imaginary. 

// a t- deformation preserves all the moments t^, k > 0, then the differential 
a{(,(){9tCdC — dtCd() extends to a holomorphic differential in D^. Take a small 
neighbourhood B of G D such that \z\ < \(\ for all 2; G B and C ^ 7) then for 
2; G B we can expand: 

dtc+(z)dz = I ( 1^ E / c'-'udc, dc] = Yl mtk)z'-'dz = 

\ k>0 "^T / k>l 

and, since (7+ is analytic in D, we conclude that 9^(7+ = 0. Therefore, the differ- 
ential C){dtCdC ~ dtCdC) is the boundary value of the differential —dtC~ {z)dz 
which has at most simple pole at 00 and holomorphic everywhere else in D'^. The 
equality 



dtto = ^ / ^(C, CMCdC - dtCdO = 

2TTt 



then implies that the residue at 00 vanishes, so dtC {z)dz is a holomorphic differ- 
ential. 

Any holomorphic differential which is purely imaginary along the boundary of a simply- 
connected domain must identically vanish in this domain. Hence we conclude that 
dtCdC — dtC,dC, — which means that the vector dtC, is tangent to the boundary and 
hence the deformation is trivial. □ 



In fact the assertion of the Proposition remains true under less restrictive conditions 
on the function a: it is enough to require that cr 7^ in some strip- like neighbourhood of 
the contour 7. 

The fact that the set of the moments is not overcomplete follows from the explicit 
construction of vector fields in the space of domains that change real or imaginary part 
of any moment keeping all the others fixed (see below). These arguments allow one to 
prove the following theorem. 

Theorem 4.1. The real parameters to, Ret^, Imt^, k > 1, form a set of local coordinates 
in the space of simply- connected plane domains with smooth boundary. 

This statement allows one to identify functionals on the space of domains D with 
functions of infinitely many independent variables to, {tk}, {ik}- 
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4.2. The Green function of the Dirichlet boundary value problem and special 
deformations. According to the Riemann mapping theorem, there exists a conformal 
bijection w : D'^ — U between D'^ and the exterior of the unit disk U = {ugC| |m|>1}. 
If the conformal map w{z) is known, one can construct the Green function of the domain 
0', 

(13) GM-iio,m^ 

which solves the Dirichlet boundary value problem in D*^. To wit, the Poisson formula 

(14) u{z) = - ((uo{OdnM^,Om 



restores a harmonic function in D'^ from its boundary value Uq = u\^. Here (9„^ denotes 
the derivative along the outward normal vector to the boundary of D with respect to the 
second variable and \d^\ is an infinitesimal element of length along the boundary. The 
Green function G{z,^) = 2^1og|z — ^| + g{z,^) in D*^ x D"^ is uniquely defined by the 
following properties [10] : 

• G{z, = G{^, z) and G{z, (!) = for any ^ e and C e 7; 

• The function g{z,^) is continuous in D'^ x D'^ and is continuous for ^ in the closure 
of D'^ at any z G D^; 

• The function g{z,C,) is harmonic in z for any G D'^ and is harmonic in ^ for any 

z e 

Let us introduce the differential operator 

(15) Viz) = do + D{z) + D{z), 
where D{z), D{z) are given by (jl]). 

Fix a point a G D"^ and consider a special infinitesimal deformation of the form 

671 

(16) San{z) = — —dn,G{a,z), ^ G 7, e 0, 

cr{z,z) 

where San{z) is the normal displacement of the boundary (positive if directed outward D) 
at the boundary point (equivalently, one may speak about the normal "velocity" of the 
boundary deformation which is Va{z) = lim^^Q{6na{z) / e) . For any sufficiently smooth 
initial boundary this deformation is well-defined as £ ^ 0. By 6a we denote the variation 
of any quantity under this deformation. 

Lemma 4.1. Let X he any functional on the space of domains D regarded as a function 
of to, {tk}, {ik}, then for any z E D'^ we have 6zX = 6'V{z)X. 

Proof. It is easy to see that 

S.to = -^^^dn^G{z,Om=e, 6ztk = -^f^^''9n,G{z,0m 

by virtue of the Poisson formula dH]). Therefore by Theorem 14.11 we have: 

dX dX dX / z~^ z~^ - \ 

'-"^ = *^*» + E a; + E ^^f' = + E — * + E — *) 

^ k>l k>l k>l k>l 

□ 
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Now we can explicitly define the deformations that change only either Xk = Hetk or 
Hk = Im tk keeping all other moments fixed. From the proof of Lemma 14.11 it follows that 
the normal displacements 6n{^) = eKe {dn^Hk{C,)) and 5n{^) = elm (c}„^iffc(^)), where 



J oo 



(the contour integral goes around infinity) change the real and imaginary parts of tk by 
±e respectively keeping all other moments unchanged. In particular, the deformation 

£71 

changes tg only. Therefore, the vector fields d/dtQ, d/dx^, d/dyk in the space of domains 
are locally well-defined and commute. Existence of such vector fields means that the 
variables tk are independent and dk = ^{dx^ —idy^,), dk = ^{^x^. + 'idy^,) can be understood 
as partial derivatives. 

We will also need the following simple lemma. 

Lemma 4.2. Let X be a functional of the form X = ^(C) C) "^(C) C) d'^C with an arbi- 
trary integrable function \E' regular on the boundary, then 

where \E'^(z) is the (unique) harmonic extension of the function \l/ from the boundary to 
the domain D'^. 

Proof. The variation of X under the special deformation fll6p is 

s^x = (f vi>(c, c)^(C, OSn.iO MCI = -BTT I ^{C, Odn,G{z, c) MCI 

Now the assertion obviously follows from Lemma [4.11 and the Poisson formula (fl^ . □ 

4.3. The dispersionless tau- function. Consider the following functional on the space 
of domains D: 

(17) ^ = ~h /{/(^(^'^)iogk"'-C"V(C'C)c^'^A. 

Theorem 4.2. ^ [121 [HI EH] j It holds 

(18) G{z, C) = ^ log \z-^ - + ^V(z)V(C)F. 

Proof. The proof consists in successive application of Lemma 14.21 and using the charac- 
teristic properties of the Green function. Applying Lemma [4.21 to f[T71) . we get: 

(19) V(^)F = -^ jj^\og\z~'-C'\<y{CX)d\ 

for z G D*^. Applying the lemma once again, we conclude that V(C)V(-2)-F is the harmonic 
continuation of the function —2 log \ z^^ — C,~^ \ from the boundary to the domain D'^. This 
function is harmonic everywhere in D'^ except ai = z, where it has the logarithmic 
singularity. It can be cancelled, without changing the boundary value, by adding the 
function 47rG(2;, (^). □ 
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Let w{z) be the conformal map from D'^ onto the exterior of the unit circle nor- 
mahzed by the conditions w{oo) = oo and w'{oo) is real positive. It has the form 

w{z) = pz + ''^^pjZ^^ , where p > 0. 

i>o 

Corollary 4.1. The conformal map w{z) is given by 



(20) 



w{z) = zexp ((-^^0 - doD{z)y^ 



Proof. From equation (fT3|) it follows that 2ttG{z, oo) = — log 1^(2;) |. Tending ^ — )• oo in 
(fT8|) and separating holomorphic and antiholomorphic parts in z, we get the result. □ 

Note that the limit z — )■ oo in ( l20l) yields logp = — | OqF. 

The following theorem establishes the embedding of the conformal dynamics into the 
dispersionless 2DTL hierarchy and identifies F with the dispersionless tau-function. 

Theorem 4.3. ( |15[ 135] ) The function F satisfies the equations 

1 _ ^-D{z)D{OF ^ ]^ ^doido+Diz)+D{0)F 

(cf m-m)- 

Proof. The proof is the same as in |15j . Combining ( !T3|) and ( !T8|) . we get 



log 



w{z) - w(0 



= log 


1 1 









iv(z)v(OF. 



1 — t(7(z)w(,^) 

Next, substituting here w{z) from (120|) and separating holomorphic and antiholomorphic 
parts in z, ^, we obtain the desired result. □ 

Important remark. Note that although the definitions of the moments and the 
function F essentially depend on the background density a, the formulas for the Green 
function and the conformal map (fT8|) . (pOj) are cr-independent. This means that the 
conformal dynamics is described by any solution to the dispersionless 2DTL hierarchy of 
this class. 

4.4. Complimentary moments. The integral of Cauchy type 

c{z) = ^ j( MiM , udC, ^ d,u{c, 0, 

defines a function which is analytic in D and D'^ with a jump across 7. Let C^{z) be the 
analytic functions defined by this integral in D and D'^ respectively. By the Sokhotski- 
Plemelj formula, the jump of the function C{z) across the contour 7 is equal to dU{z,z): 

(21) iC^z)-C-iz))\^^^ = dUiz,z). 



As it follows from the proof of Proposition 14. ![ C^{z) is the generating function of the 
moments tk'- 
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Similarly, C {z) is the generating function of the set of complimentary moments Vk'- 

•'^ ^ k>i 
which are given by the integrals 

Vk = — (f z^dU {z, z)dz = - [ [ z^a{z, z) £z , k> 1. 

They are functions of the moments to, {tk}, {ik}- It is also useful to introduce the loga- 
rithmic moment 

vq = — / / \og\z\'^a{z,z) (fz. 



71" „ „ ^ 

Theorem 4.4. (cf. [HI [331 [35] j The following relations hold: 

(22) vo = doF, Vk = dkF, Vk = BkF, k>l. 

Proof. Applying Lemma [4.21 to ([TT]) . we get ( !T9|) . The expansion of both sides in powers 
of 2, z yields ([22]). □ 

Proposition 4.2. Suppose that only a finite number of the moments tk are different from 
0. Then the tau-function / f77|) can be represented as 



(23) 2F = -^ll Uad^z + t^v^ + ^{tkVk + t 



tkVk}- 

k>l 



Proof. Set 0(z, z) = — ^ fj^log \z ^ — C ^ \ cr(C) C) d'^C^ then under the assumption of the 
proposition we have the expansion 

[z,z) = -U{z,z) +to\og\z\'^ + '^{tkz'' + ikz'') 



k>l 



valid everywhere in D \ {0}. Substituting this into f[T7j) written in the form 
2F = i JJ^(j)ad'^z, performing the termwise integration and using the definition of the 
complimentary moments, we get ([23]) . □ 



5. CONFORMAL DYNAMICS IN A SYMMETRIC BACKGROUND 

5.1. The general axially symmetric case. The case when the background density 
function a (and the function U) is axially symmetric, i.e., depends only on \z\, is of a 
special interest. Note that for any axially symmetric background it holds zdU = zdU. 
We denote the symmetric density by (T(|2;p). In this case the conformal dynamics provides 
a symmetric solution to the dispersionless 2DTL, so the function F has the explicit Taylor 
expansion obtained in section 2. Indeed, the symmetry implies that when all moments tk 
at A; > 1 are equal to 0, the domain D is just a disk of radius R such that 

(24) to = a{x)dx, vq = / log xa{x)dx 

Jo Jo 

and the complimentary moments Vk with k > 1 are zero. From (122 p we see that 

(25) Vkl^ = dt.Fl^ = forallfc>l 
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so the solution is indeed symmetric in the sense of section 2. One can see that 
dvo/dtf) = logi?^, i.e., p = 1/R as it should be. We also have 

rR'^ ry 

(26) F\^^ = a{y)dy \ogxa{x) dx, /(to) = e^l"o = R'ito), 

where i?^ should be understood as a function of to implicitly given by (124|) . 
The Green function at = is 



271^(^,0= log 



R 



log|.-^-r^|+logi?-Re5^ — 

k>l ^ 



zC-R^ 
Comparison with f fTSj) yields 

(27) d,^d,,Fl^ = 0, d,^\Fl^ = kR'%, 

in accordance with Lemma [2.11 

Note that in the axially symmetric case the moment Vq can be represented as a contour 
integral: 



— (f {\og\z\'^ dU - z-^U) dz. 



2 

5.2. Examples of symmetric solutions. 

5.2.1. The homogeneous density. An important example is the homogeneous density 
(28) a{z, z) = {z-zr-\ U{z, z) = ^ {zzT 



a 



with some a G M. Assume that a > 0. In this case to = -R^°/a, so /(to) = (ato)^^° and 
(recall that for symmetric density f = R 



2^ 



Proposition 5.1. Assuming that only a finite number of the moments t^ are different 
from 0, the dispersionless tau-function for this solution is quasi-homogeneous, that is it 
obeys the relation 

(29) 4aF = -tl + 2ato9oF + ^(2« - k) {tkdkF + iAF) . 

k>l 

Proof. We use equation (123|) . In the integral term we write Ua = d{UdU) — dU dU and 
notice that for the particular function U we have Ua = dUdU, and also zdU = aU, so 
Ua = I d{UdU). This allows us to transform the 2D integral to a contour integral: 

- // Uad^z=^ LzdUf — . 
-KjJo A-KiaJ^ z 

Now recall ( 12T|) and represent dU = — (on 7), with being a polynomial. 
Shrinking the integration contour to 00, we obtain: 



-1. fizdUr- = tl + 2j2khVk. 



k>l 

Since the initial integral is obviously real, we conclude that X]fe>i ktkVk is a real quantity. 
The quasi-homogeneity relation fl2^ follows. □ 
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Remark. The case a = 1 {a{z,z) = 1) was considered in [T8l 133] . The Taylor 
coefficients for this solution were found in [22j. In this case the conformal dynamics is 
identical to the Darcy law for the motion of interface between viscous and non-viscous 
fluids confined in the radial Hele-Shaw cell, assuming that there is no surface tension at 
the interface. The normal velocity of the interface 7 is proportional to the gradient of 
the Green function of the Laplace operator. This type of processes is known as Laplacian 
growth (see, e.g., [TJ [2Q]). As is mentioned in |21] . the conformal dynamics in the case 
a = 1/N for an integer N > can be mapped to a Laplacian growth process in a sector 
with angle 2it/N and periodic conditions at the boundary rays (a cone). With proper 
modifications, the case a < (and, in particular, a = —1) can be also considered and 
mapped to the Laplacian growth in the compact interior domain bounded by the curve 

7- 

5.2.2. The solution yelding the Hurwitz numbers. In this case 

1 1 r zz^'^ 

(30) a(z, z) = ^, Uiz, z) = — log - 

which can be formally regarded as a limit a — )■ of f l28p but the limit is tricky and this 
case deserves a separate consideration. In this case 

i^L„ = ^ + f logQ, m = QeP^^. 

Remark. As is shown in [34j, the conformal dynamics in this case can be mapped to a 
Laplacian growth process in an infinite channel with periodic conditions in the transverse 
direction, i.e., on the surface of an infinite cylinder of radius 

Proposition 5.2. (]34J j The dispersionless tau-function for this solution obeys the fol- 
lowing homogeneity relation: 

2F = TdrF + t^doF + {tkdkF + tAF) , r = . 

k>l 

1 f / 1 ^ I ^ \ dz 1 f / \z\'^ \ dz 

Proof We have h = log^— j -k>l.t, = ^^f^ log(— j -, so that 

the variables tk = (3tk are /3- independent. As is seen from equation (fTTl) with a as in f l30|) . 
F is of the form F = P~'^F, where F is /3-independent. Therefore, 

F(/3, {t,}) = F(/3, {4//3}) = p-'F = p-'F{l, {4}). 
Taking the total /3-derivative of the identity (3'^F{I3, {4//?}) = F(l, {4}), we get: 
-2/3^ F - I3%F + (to doF + 2 Re^ 4 -9^^) = 

k>l 

or 

2F = -l3dpF + to doF + J] (t^ d^F + h BkF^ . 

k>l 

□ 



The Taylor expansion of this function is given by (fT2|) . with the Taylor coefficients being 
esseentially the double Hurwitz numbers for connected genus coverings. The homogene- 
ity property can be explicitly seen from the Taylor expansion. 
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Proposition 5.3. The dispersionless tau-function for the solution determined by the data 
satisfies the relations 

OF tl 



'^0 



dlogQ 2 
(31) 

dF t^ 1 

Qo = + toY.ktkdkF + [kltkk dk+iF + {k + l)tk+i duF diF) . 

^ k>l k,l>l 

Proof. The first formula is obvious from the structure of the Taylor expansion f fT2|) . An- 
other proof, similar to that of Proposition 15. ![ is given in [3^. For the proof of the second 
formula we note that in the case (130|) f/cr = | d{UdU) and thus 

Proceeding as in the proof of Proposition I5.H we obtain 

j) {zdUf — = -| + ktk ^fc + 2 5Z {k^^kkVk+i + {k + l)tk+iVkVi) . 

"^T k>l k,l>l 

Then from Proposition 14.21 we have: 

2F = toVo + ^itkVk + tkVk) ^ - /3to^ ktkVk 

k>l k>l 

(32) 

_ ^ 
2 



[kltktiVk+i + {k + l)tkMVkVi 



kl>l 



Combining this with the homogeneity property (Proposition 15. 2p . we arrive at the second 
formula in (IHT]) . □ 

Remark. The double sum in the second equation in (13T]) is the genus part of the 
celebrated cut-and-join operator [5j, see also [27J. 
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